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Abstract 

Local magnetic relaxation measurements in YBa2Cu307_ x crystals show evi- 
dence for plastic vortex-creep associated with the motion of dislocations in the 
vortex lattice. This creep mechanism governs the vortex dynamics in a wide 
range of temperatures and fields below the melting line and above the field 
corresponding to the peak in the "fishtail" magnetization. In this range the 
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activation energy U p i, which decreases with field, drops below the elastic (col- 
lective) creep activation energy, U e i, which increases with field. A crossover 
in flux dynamics from elastic to plastic creep is shown to be the origin of the 
fishtail in YBa2Cu307_ :r . 

PACS numbers: 74.60.Ge, 74.72.Bk 
Magnetic relaxation in high-temperature superconductors has been traditionally de- 
scribed in terms of collective vortex-creep based on the concept of elastic motion of the 
vortex lattice This approach successfully explained wide range of vortex dynamics 

phenomena. On the other hand, under certain conditions effects of vortex lattice plasticity 
may become the dominant factor that determines the vortex dynamics ||. For example, 
numerical simulations of strongly pinned vortex system reveal vortex motion dominated by 
plastic deformations ||. Experimentally, effects of plastic vortex behavior were observed 
in transport measurements in close vicinity of the melting transition where the pre-melting 
softening of the lattice enhances the role of spatial inhomogeneities, resulting in the tearing 
of the vortex lattice under applied currents |5|J§]. In this paper we demonstrate that plastic 
deformations dominate vortex-lattice motion in the creep regime in YBa2Cu307_ ;c (YBCO) 
crystals far below the melting transition, in the region that was believed to be governed by 
elastic motion. We find that this plastic motion governs the flux creep in YBCO crystals 
at elevated temperatures at fields above the characteristic field B p corresponding to the 
peak in the 'fishtail' magnetization [pHbj; and it affects the corresponding shape of the 
magnetization curves. 

The plastic vortex motion can be classified into three main categories, i) Vortex channel- 
ing along easy paths in the pinning relief in between rather stationary vortex-lattice islands. 
Such behavior was observed in numerical simulations in presence of very strong pinning [|J . 
ii) Vortex motion that resembles ice floe in which large pieces of vortex-lattice slide with 
respect to each other. As was simulated numerically these two modes of dynamics are 



believed to be the cause of the highly unstable character of resistivity in the close vicinity of 
the melting transition [f||§. iii) Dislocation mediated plastic creep of vortex lattice similar 



to diffusion of dislocations in atomic solids ||16|| . This type of vortex behavior was consid 



ered previously |IT]IB| without an experimental evidence. The observations described below 
strongly suggest that this mode of plastic vortex creep governs magnetic relaxation over a 
substantial part of the YBCO magnetic phase diagram. 

Local magnetic relaxation measurements were performed on a 1.2 x 0.5 x 0.3 mm 3 single 
crystal of YBCO (T c ~ 91 K) using an array of microscopic GaAs/AlGaAs Hall sensors with 
30 x 30 fim 2 active area and sensitivity better than 0.1 G. The probes detect the component 
B z of the field normal to the surface of the crystal. Temperature stability and resolution 
was better then 0.01 K. 

After zero-field-cooling (zfc) the sample from above T c to the measurement temperature 
T we measured the full hysteresis loops for all the probes. The first field for full penetration 
H* was measured directly by the probe at the center of the sample. After repeating the 
zfc process, a dc field H was applied parallel to the c-axis and the local induction B z was 
measured at different locations as a function of the time t for an hour. These relaxation 
measurements were repeated after the field was increased by a step AH > 2H* up to the 
irreversibility field H irr or the maximum field of the experiment (1.6 T). 

The inset to Figure 1 shows typical hysteresis loops, B|' — H vs. H&tT = 85K for four 
probes (i = 5, 6, 7, 8 located at 70, 130, 190, and 250 /xm from the edge towards the center). 
Each probe exhibits a clear fishtail behavior with a maximum in local magnetization at field 
Bp ~ 0.4 T. The width of the loop is largest in the center of the crystal and decreases 
towards the edges, as expected from basic considerations based on a modified Bean model 

In Figure 1 we show the time evolution of the gradient (Bf^ - Bi 7 A /Ax between h = 8 
sec and t 2 = 3600 sec, as a function of the applied field. In our geometry (aspect ratio 
= 3/5) this gradient is proportional to the persistent current density J that can be readily 
evaluated using sensors 6 and 7 which are located not too close to the center or the edge 
of the crystal [fZ0}-|2"2"|] . Note that the position of the fishtail peak B p shifts from 0.42 T 



to 0.34 T during the relaxation. The total relaxation of the current, AJ, during the time 



window of the measurement is large for B < B p and it increases to even larger values above 
B p . Furthermore, the relative change of the persistent current, AJ/J, is also large (e.g. 
A J / J pd 0.8 at H = 0.9 T), implying that J <C J c on both sides of the fishtail peak. These 
strong relaxations imply that dynamic effects determine the shape of J(B) and rule out 
the possibility [pj,|5||n2"| that either branch of the fishtail is determined by the critical current 
density J C (B). 

Knowledge of the time and spatial field-derivatives, dB z /dt and dB z /dx, enable direct 



determination [2^] of the activation energy U(J, B) associated with the flux creep by using 
the diffusion equation: dB z /dt = —d/dx(B z v), where the effective vortex velocity v is 
proportional to exp(—U/kT). Typical U vs. J data, at 85 K, are shown in Fig. 2 for 
different fields, in the range of 0.05 T to 0.8 T. This figure exhibits a dramatic crossover in 
the slope \dU/dJ\ around B p = 0.4 T. In order to quantify this crossover, we start by using 
the prediction of the collective creep theory for J <^ J c |I|,|17| : 

U(B, J) = U (B)(J C /Jr oc B v r\ (1) 

where the positive critical exponents v and /i depend on the specific pinning regime. We 
note that the range of the experimentally accessible ?7-values in Fig. 2 is almost independent 
of the field p3| . This implies J oc B u ^, i.e. J grows with field for J J c - Obviously, the 
collective creep dynamics cannot explain the decrease of J with B observed above B p . 

The inset to Fig. 2 shows the /i-values obtained by fitting the U(J) data, using Eq. [p. 
At low fields /i ~ 1 and it then increases to /i ~ 2 just below B p . This confirms || that 
below the peak the relaxation is well described by the collective creep theory. The latter 
predicts /x = 1 in the intermediate bundle regime and fi = 5/2 in the small bundle regime 
|24| . However, above the peak, [i drops sharply to values below 0.2. Within the collective 



creep theory this would imply an inconceivable crossover to a single vortex regime {\x = 1/7) 
which is expected only for low fields and high values of J Thus the /i-values above B p 
are inconsistent with the collective creep theory. 

The failure of the collective creep theory in explaining the data above B p can be further 



demonstrated by analyzing the exponent v in Eq. ([I]). Separation of variables in this equation 
implies that a smooth U(J) function can be obtained by proper field-scaling of the data sets 
at various fields. However, as demonstrated in Fig. 3 we need two exponents of opposite 
signs, v = 0.7 and v = —1.2, to scale the data of Fig. 2 for fields below and above B p , 
respectively. In the collective creep theory U increases with B, thus the negative ^-value is 
inconsistent with this theory. 

We now show that the experimental data above B p can be well explained by a plastic 
creep model based on dislocation mediated motion of vortices similar to diffusion of disloca- 
tions in atomic solids [ 16| . The activation energy U pl at J = for the motion of a dislocation 



in the vortex lattice can be estimated as 18 



U pl (B) ~ ee a oc 1/Vb (2) 



where e is the vortex line tension, e = \/m a b/m c is the anisotropy parameter, and 



•)/B is the mean intervortex distance. This estimation assumes a formation of a 
dislocation semi-loop between two valleys separated by a distance a fl6|]. One notices that 



U pl decreases with field in contrast to the collective creep activation energies U e i, which 
increases with field, see Eq. (||). Of course, the creep process is governed by the smaller be- 
tween U e i and U p i. Thus, at low fields where U p i > U e i, the latter controls the flux dynamics. 
But, as B increases and U p i becomes less than U e i, a crossover to the plastic creep regime is 
expected. 

The values of U pl (B) can be extracted from the measured U(J) curves of Fig. 2 assuming 
an expression for U p i(J) taken from the dislocation theory [|16[], substituting the current 
density for the strain: 

U pl (J) = U° pl (l - \fJJjf) , (3) 

where Jf is the critical current which corresponds to the plastic motion. The derived 
values are shown in the inset of Fig. 3 as a function of B. The solid line in the inset is a 
power-law fit to the experimental data, U pl cx B~ 7 . Clearly, U pl decreases with the field, 
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although with an exponent —0.7 rather than the expected —0.5. As we show below, the 
same exponent (—0.7) is also involved in determining the temperature dependence of B p . 

The fishtail peak location B p can be determined from the condition U e \ = U p i for the 
same J. Note that U p i = U pl since J <C Jf. Using the logarithmic solution of the flux 
diffusion equation |T[ U e i = kThi(t/t ) together with Eq. (|2|), we get B p oc l/ln 2 (t/t ), i.e. 
the peak position should shift with time towards low fields, as observed in Fig. 1. Similarly, 
for the temperature dependence of B p one obtains: 



B p oc el oc 1/A 4 oc (l - (T/T c ) 4 ) . (4) 

A fit of this expression to the experimental data yields a modest agreement. However, 
as noted above, in fact U pl oc B~ 07 , thus B P (T) oc e^ ' 7 , i.e. B p oc (1 — (T/T C ) 4 ) 1A . 
Indeed, a perfect fit (B p , solid line in Fig. 4) is obtained with this expression. Note that 



a similar variation of the exponent was found in some experiments [25] for the melting 



line B m (T) oc (1 — (T/T C )) L4 , while the theory |T[ predicted a critical exponent of 2. This 
interesting similarity may support previous claims |J that the plastic motion of defects in 
the vortex lattice is a precursor to the melting transition. 

In Figure 4 we show that the plastic creep regime in the B — T phase diagram (shaded 
area) covers a substantial area between the B P (T) line found in this experiment and the 



melting line B m {T) described in Ref. |25 



In conclusion, our data clearly indicate two different flux-creep mechanisms above and 
below the peak in the magnetization curves. In particular, flux creep with activation energy 
decreasing with field plays an important role above the peak. The data in this regime cannot 
be explained in terms of the traditional collective creep theory based on the concept of elastic 
motion of the vortex lattice. On the other hand, the data show good agreement with the 
dislocation mediated mechanism of plastic creep analogous to plasticity in atomic solids. 
This observation leads also to the conclusion that the origin of the fishtail in YBCO crystals 
is a crossover from elastic to plastic creep. The predictions of this model for the time and 
temperature dependence of the location of the peak are well confirmed in the experiments. 
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FIGURE CAPTIONS 

Fig. 1 Relaxation of persistent current J calculated from the field gradients for different 
applied fields at T = 85 K. Note the shift in the peak position B p during the exper- 
imental time window 8 — 3600 sec as noted by the arrows. Also note the increase in 
the relaxation rate above B p . Inset: Local hysteresis loops for four Hall probes vs. 
applied field H at T = 85 K. The width of the loops increases from the edge of the 
sample towards the center. 

Fig. 2 U vs. J for the indicated fields (0.05 T — 0.8 T) at T = 85 K. The lines are guide to the 
eye. Note the change in the slope dll/dj above and below the peak field B p ~ 0.4 T. 
Inset: The critical exponent /x, see Eq. (|1|), as a function of field, at T = 85 K. 

Fig. 3 Scaling of U(J, S)-curves below and above B p atT = 85 K. Above the peak U oc B 7 
and below the peak U oc B~ 12 . The inset shows the derived U pl vs. B. The solid line 
is a fit to U p i oc B~ 7 . 

Fig. 4 Vortex-creep phase diagram for YBCO. The plastic creep regime is limited between 
B p (T) (solid line) and the melting line B m (T) (dashed line, taken from [^J). Data 
points in the B p (T) curve were determined in this experiment at t\ = 8 sec. The solid 
line is a fit to B p oc (1 — (T/T C ) 4 ) L4 . 
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Fig. 2. Abulafia et al. J [kA/CHl ] 




5000 1 10 4 1.5 10 4 

Fig. 3. Abulafia et al. J [ kA/dll ] 



